Quantum phases of dipolar bosons in optical lattices 
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The ground state of dipolar bosons placed in an optical lattice is analyzed. We show that the 
modification of experimentally accessible parameters can lead to the realization and control of dif- 
ferent quantum phases, including superfluid, supersolid, Mott insulator, checkerboard, and collapse 
phases. 
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The Bose-Einstein condensation (BEC) of dilute 
atomic gases jl| has opened a new interdisciplinary area 
of modern atomic, molecular and optical (AMO) physics 
on one side and condensed matter physics on the other: 
the study of ultracold weakly interacting trapped quan- 
tum gases 0. So far most of the experiments in this 
area have been very accurately described by the semi- 
classical mean-field method and its extensions, based 
on the Gross-Pitaevskii (GP) and Bogoliubov-de Gcnncs 
equations J3|. However, experimental techniques have 
recently progressed to a stage at which mean-field meth- 
ods cease to provide an appropriate physical description. 
In this sense, experiments on Feshbach resonances at 
JILA Q allow the modification of the s-wave scatter- 
ing length to such large values that the mean-field pic- 
ture is no more applicable. Similarly, the achievement 
of BEC in metastable helium Q opens the possibility to 
study higher-order correlation functions, whose analysis 
requires theoretical approaches beyond mean field. The 
recent observation of the Mott insulator-superfluid phase 
transition in ultracold atomic samples in optical lattices 
P|, predicted in 0, belongs to the same category, but 
at the same time initiates a new research area of AMO 
physics: the physics of strongly correlated quantum gases. 
The experiments of || are relatively easy to accurately 
control and manipulate and thus provide a novel and par- 
ticularly promising test ground for theories of quantum 
phase transitions |S|] , which have traditionally dealt with 
condensed-matter systems rather than with atomic gases. 

The influence of dipole-dipole forces on the properties 
of BEC has also drawn a considerable attention recently. 
It has been shown that these forces significantly mod- 
ify the ground state and collective excitations of trapped 
condensates |)|-|lT|. Dipole-dipole interactions are also 
responsible for spontaneous polarization and spin waves 
in spinor condensates in optical lattices JL^] and may lead 
to self-bound structures in the field of a traveling wave 
p3| . In addition, since dipole-dipole interactions can be 
quite strong relative to the short-range (contact) inter- 
actions, dipolar particles are considered to be promis- 
ing candidates for the implementation of fast and robust 
quantum-computing schemes |14[|15[|. Sources of cold 



dipolar bosons include atoms |16| or molecules [n7[ with 



permanent magnetic or electric dipolc moments. Other 
possible candidates could be atoms with electric dipoles, 
induced either by large dc electric fields || or by opti- 
cally admixing the permanent dipole moment of a low- 
lying Rydberg state to the atomic ground state in the 
presence of a moderate dc electric field [[ll 14 . 



This Letter is devoted to the analysis of the ground 
state of an ultracold gas of polarized dipolar bosons in 
an optical lattice. The ground state of a gas of short- 
range repulsively-interacting bosons in a periodic poten- 
tial can be either in a superfluid phase or in a Mott- 
insulating phase, characterized by integer boson densities 
and the existence of a gap for particle-hole excitations 
The superfluid-Mott insulator transition in cold 
bosonic atoms in optical lattices has been recently the- 
oretically analyzed and experimentally demonstrated 
H . For the case of finite-range interactions new quan- 
tum phases have been predicted |l9[ ] , including supersolid 
phases which combine both diagonal and off-diagonal 
long-range ordering. To the best of our knowledge, 
dipole-dipole interactions have not yet been discussed 
in this context. We show in the following that these 
interactions, which are long-range and anisotropic, lead 
to new interesting properties. The long-range character 
of the dipole-dipole potential provides a rich variety of 
quantum phases. Moreover, we show that the interac- 
tions in a gas of dipolar bosons are easily tunable, allow- 
ing for the experimental engineering of quantum phase 
transitions between various kinds of ground states. Such 
a highly controllable system may be crucial in answer- 
ing some unresolved questions in the theory of quantum 
phase transitions (e.g. the existence of a yet-unobserved 
supersolid |^0| , or a Bose metal at zero temperature |2lJ ) . 

A dilute gas of bosons in a periodic potential (e.g. in 
an optical lattice) can be described with the help of the 
Bose-Hubbard (BH) model 0. For particles interacting 
via long-range forces the BH Hamiltonian becomes: 
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where hi is the annihilation operator of a particle at the 
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lattice site i, which is considered as being in a state de- 
scribed by the Wannier function w(r — r*) of the lowest 
energy band, localized on this site. This implies the as- 
sumption that the energies involved in the system are 
small compared to the excitation energies to the second 
band. We denote the position of the local minimum of the 
optical potential as r,;, and the number operator for the 
site i as rij = b\bi. In Eq. (|l|) only the nearest-neighbor 
tunneling is considered, which is described by 

J = [ w*(r - ri )[-|-V 2 + V,(r)Mr - rj ) d 3 r , (2) 
J 2m 

where j and i are indices of the neighboring sites, and 
^( r ) = Z~2^= x y z cos 2 (fc^^) is the optical lattice po- 
tential with the wavevector k. The interparticle interac- 
tions are characterized by the parameters 

U a = J |«>(r-r 4 )|V int (r-r')Kr'-r i )| 2 d 8 r dV , 

(3) 

where |r^ — Tj\ = 47rcr/|k|. Uq determines the on-site 
interactions, U ai the nearest-neighbor interactions, U a2 
the interactions between the next-nearest neighbors, etc. 
Consequently, the respective summations in Eq. (Q) must 
be carried out over appropriate pairs of sites which are 
marked by < . . . > for the nearest neighbors, << ... >> 
for the next-nearest neighbors, etc. In the 2D calcula- 
tions presented below, we have taken into account inter- 
actions with up to 4 neighbors (a% = 1, = V2, 03 = 2, 
04 = v5), since in the particular cases we analyzed the 
effects of interactions of a longer range are negligible. In 
the case of polarized dipoles the interaction potential is 
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where the first term is the dipole-dipole interaction char- 
acterized by the dipole d and the angle 9 between the 
dipole direction and the vector r — r', and the second 
term is the short-range interaction given by the s-wave 
scattering length a and the atomic mass m. 

We find the ground state of the system using a varia- 
tional approach (see (7) and references therein) based on 
the Gutzwiller ansatz \^mf) = Ili I&) f° r the ground- 
state wavefunction, where the product is over all lattice 
sites. The wavefunctions \<pi) for each site are expressed 
in the basis of Fock states, \<j>i) = Y^=o fn\ n )i> where n 
indicates the occupation number. The coefficients {/^} 
are found by minimizing the expectation value of the 
Hamiltonian (|l]) in the state \^mf) under the constraint 
of a fixed chemical potential fi. 

In the following we consider ID and 2D geometries. 
Low-dimensional BECs have been achieved in recent ex- 



periments j 22 by transversally confining a condensate in 
a tight optical or magnetic harmonic trap. A ID or 2D 



lattice is created by a laser standing wave, which gener- 
ates a periodic optical potential (^3|. We have carried 
out the minimization of (^ mf\H — fiJ2 i n^MF) for ID 
lattices with up to 20 sites and for square 2D lattices 
with up to 9 x 9 sites. A similar qualitative picture is 
expected for a larger number of lattice sites. Since in 
ID and for systems with few atoms the application of a 
mean-field calculation could be questionable (due to the 
possibly important role of fluctuations), we restrict our 
discussion of the BH Hamiltonian (|l|) to the 2D case. 

For a square 2D lattice in the xy plane, the wave- 
functions \(f>i) can be written as a product of Wannier 
functions in the x and y directions and Gaussian func- 
tions in the z direction. There are two generic situations 
for dipoles in 2D lattices, namely (i) the dipole is along 
the z direction or (ii) the dipole direction is in the xy 
plane. This follows from the fact that two dipoles expe- 
rience maximal attraction along the dipole direction and 
maximal repulsion in the transversal plane. As shown 
in Ref . |ll]] , the mean-field dipole-dipole energy critically 
depends on the shape of the bosonic cloud, which can 
be altered by modifying the trap. It is intuitively clear 
that a cloud elongated in the dipole direction is unsta- 
ble due to the predominance of attractive interactions. 
On the contrary, the cloud may be stable if it is broader 
in the transversal plane than in the dipole direction. In 
particular, for spherically symmetric wavefunctions \4>i) 
the on-site averaged dipole-dipole potential vanishes, and 
only the short-range interactions contribute to Uq. 

Let us first focus on the case (i). Since the dipole- 
dipole mean-field critically depends on the shape of the 
cloud, the balance between attractive and repulsive in- 
teractions can be easily manipulated either by modifying 
the wavelength and intensity of the lattice or by chang- 
ing the transversal trapping. In the following we employ 
the latter possibility to provide an example of how dif- 
ferent phases of the BH Hamiltonian ([j]) may be accom- 
plished just by changing the magnitude of controllable 
external fields. The expectation values (bi) provide the 
superfluid order parameter. It is non zero and constant 
for all lattice sites in the superfluid phase, whereas it is 
periodically modulated in the supersolid phase. Fig. |l| 
shows the maximal (circles) and minimal (squares) value 
of I (bi) I and of the occupation number (m) as a function 
of the aspect ratio of the on-site wavefunction \<fii). The 
aspect ratio is defined as the square of the ratio between 
the width of \<fii) in the x direction and the width in the z 
direction, L = (l x /l z ) 2 . In the calculations presented in 
Fig. |l|, we have considered the case of 23 Na atoms with 
an induced dipole moment of 0.334 Debye, placed in a 
lattice of wavelength A = 795 nm. The maximum of the 
lattice potential is V® y — 10E r , where E r = h 2 k 2 /2m is 
the recoil energy. In our simulations we fix the chemical 
potential fi — 0.082 E r , which determines the mean num- 
ber of atoms and mean density. The quantum phases ap- 
pearing for the different aspect ratios depend, of course, 
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on the chosen physical parameters. However, we have 
observed a similar picture of tunable quantum phases for 
every set of parameters that we have considered. 
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FIG. 1. Maximal (circles) and minimal (squares) values of 
the superfluid parameter (bi) and of the occupation number 
(rii) as a function of the aspect ratio L of the on-site wave- 
functions. 

For L > 1.1 a checkerboard insulating phase is 
achieved, in which a site occupied by exactly one atom 
is followed by an empty site. This phase is a result of 
the long-range repulsion between particles in the pres- 
ence of a relatively weak tunneling, which prevents the 
appearance of a superfluid. For 0.9 < L < 1.1 the sys- 
tem enters a supersolid phase, possessing both diagonal 
and off-diagonal long-range order, in which the system is 
superfluid, but the superfluid parameter shows a slight 
periodic modulation. For 0.7 < L < 0.9 the influence of 
tunneling relative to the long-range interactions is large 
enough to enforce a homogeneous superfluid phase. For 
0.57 < L < 0.7 the system is a supersolid with a strongly 
modulated superfluid parameter. This phase appears due 
to a significant mutual cancellation of the on-site dipole- 
dipole interactions (attractive for L < 1) and the always 
repulsive short-range potential - the system enters an 
interesting purely long-range regime with the local inter- 
actions essentially absent. In such a case the ratio J/Uq, 
which governs the insulator-superfluid crossover p8[ , in- 
creases, driving the system from an insulating phase to 
a superfluid one. On the other hand, the long-range in- 
teractions, characterized by the coefficients U ai , remain 
considerably large and positive. As a consequence, a pe- 
riodic modulation of the superfluid parameter occurs. Fi- 



nally, for L < 0.57 the system undergoes local collapses 
due to the attractive local interactions. This last regime 
has been confirmed by checking that the maximal possi- 
ble occupation per site is always achieved, independently 
of the value of such maximal occupation. 
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FIG. 2. Maximal (squares) and minimal (circles) values of 
the superfluid parameter (6j) and of the occupation number 
(rii) as a function of the tilt angle a. 

Another simple experimental control knob is provided 
by the angle a between the dipole direction and the vec- 
tor normal to the 2D lattice plane. For a = we recov- 
ered the previous results. For a > the coefficients U a 
depend not only on the distance between neighbors but 
also on the angle between the projection of the dipole 
direction on the lattice plane and the vector joining the 
corresponding sites. In Fig. || we present a sequence of 
quantum phases obtained when the angle a is varied. For 
this calculation, the aspect ratio is fixed to L = 0.5 and 
the rest of the parameters is kept the same as those of 
Fig. |l|. For a approaching ir/2, we observe only Mott 
insulator phase or superfluid one, since when the projec- 
tion of the dipole onto the lattice plane is sufficiently large 
the dipole-dipole on-site interaction becomes positive and 
reinforces the repulsive on-site contact interactions. In 
other words, tilting the dipoles towards the lattice plane 
brings the system back to a situation of dominant local 
interactions J7|JT^]. For a approaching zero, we observe 
collapse in this particular case, as expected for L = 0.5 
from Fig. |l|. Additionally, we should point out that we do 
not observe anisotropic phases, which one might expect 
due to the anisotropy of interactions, since we employ 
periodic boundary conditions in our simulations. 
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Finally, let us stress that in the analysis of the BH 
Hamiltonian ([!]) we have expanded the field operator in 
the basis of Wannier functions, which are exact solutions 
of the single-particle problem in a periodic potential. 
This method should give correct results as long as the 
mean occupation of sites is of the order of unity. However, 
nowadays it is possible to load large Bose-Einstein con- 
densates into optical lattices |^3| , which results in a very 
high occupation of sites. In such a situation, as long as 
the lattice potential does not prevent the establishment 
of a common phase between sites, the GP equation, rou- 
tinely used to describe the condensate wavefunction in 
harmonic traps, should provide a correct description of 
the system Jll],|24|]. In the presence of sufficiently strong 
dipole-dipole forces we can neglect the short-range inter- 
actions and the time-independent GP equation reads: 

h 2 

^{r) = {-—V 2 + V l (v) + V t (r) 

+ d 2 |^ 1 k 3C r ° / ^ |^(rO| 2 Wr), (5) 

where ip( r ) is the wavefunction of the condensate (nor- 
malized to the total number of particles N) and Vt(r) = 
mu> 2 r 2 Jl is a spherically symmetric harmonic trap with 
frequency ui. In the absence of the lattice potential U(r), 
the condensate is stable as long as a = N^^y^^d 2 does 
not exceed some critical value a cr (for a spherical trap 
o'er = °"cr ~ 4.3 pH ). We have observed that, by raising 
various ID and 2D lattice configurations, one can either 
destabilize the condensate for a < a® r or make it stable 
in the regime a > a® r . For instance, the BEC is stabi- 
lized for a ID lattice whose wavevector is along the dipole 
direction, or a 2D lattice on a plane which contains the 
dipole direction. The discussion of these results will be 
presented in detail elsewhere. Let us just stress here that 
dipolar gases provide also in this situation a unique and 
very efficient possibility of coherent control of a BEC. 

In this Letter we have analyzed the ground state of 
dipolar bosons placed in an optical lattice. We have 
shown that by modifying well-controllable parameters, 
different quantum phases can be accomplished, includ- 
ing superfluid, supersolid, Mott insulator, checkerboard 
and collapse phases. 
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